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A Priori Error Estimates for Mixed Einite Element d-Schemes for the 

Wave Equation* 

Samir Karaa^ 

Abstract 

A family of implicit-in-time mixed finite element schemes is presented for the numerical approxi¬ 
mation of the acoustic wave equation. The mixed space discretization is based on the displacement 
form of the wave equation and the time-stepping method employs a three-level one-parameter 
scheme. A rigorous stability analysis is presented based on energy estimation and sharp stability 
results are obtained. A convergence analysis is carried out and optimal a priori L°°{L^) error 
estimates for both displacement and pressure are derived. 

Key words - wave equation, mixed finite elements, error estimation, energy technique 
Mathematical subject codes 65L05, 65M12, 65M60, 65M15 

1 Introduction 

The acoustic wave equation is used to model the effects of wave propagation in heterogeneous media. 
Solving this equation efficiently is of fundamental importance in many real-life problems. In geophysics, 
it helps for instance in the interpretation of the seismograph field data and to predict damage patterns 
due to earthquakes. Using finite element methods for its approximation is attractive because of the 
ability to handle complex discretizations and design adaptive grid refinement strategies based on error 
indicators. 

Previous attempts on wave simulation by finite elements have used continuous Galerkin methods 
pm El 01181122], mixed finite element methods [61 mmiisiisni 1111126], and discontinuous Galerkin 
methods uni n m [25]. In a mixed finite element formulation both displacements and stresses are 
approximated simultaneously. This approach provides higher-order approximations to the stresses. 
This property is important in many problems, in particular in modeling boundary controlability of 
the wave equation m- One of the main difficulties of the mixed finite element techniques is the 


requirement of compatibility of the approximating spaces for convergence and stability. 

Given a bounded convex polygonal domain 12 in R”^, m = 2, 3, with boundary dQ = T/) U T^r, 
and unit outward normal v, the general form of the wave equation is 

putt + V-f = f inf2x(0,r), (1.1) 

V-tt = 0 onr£,x(0,T), (1.2) 

uu = 0 onr7vx(0,T), (1-3) 

u(-,0) = inn, (1.4) 

Ut{-,0) = inn, (1.5) 
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where u is the displacement, p is the density, and r is the stress tensor given by the generalized Hooke’s 
law T = A(V • u)I + /i(Vit + (Vw)^). Here A > 0 and p are the Lame coefficients characterizing the 
material. The function / represents a general source term and vP and are initial conditions on 
displacements and velocities. We assume that /, vP and are smooth enough so that there is a 
unique solution u € C^((0,T) x fl) to (ll.ll) - (ll.5|) , see m- 

The limiting case of (inD with /i = 0 is referred to as the acoustic wave equation, which is 

putt + y ■ {X{V-u)!) = f. (1.6) 

It is assumed that p and A are bounded below and above by the positive constants po; Pi-, and 
Ai, respectively. This vector equation is equivalent to the scalar wave equation after making the 
substitution p = AV • u. The mixed method is established by using this relationship, leading to the 
coupled system 


pUtt-^P = f inflx(0,r), (1.7) 

\~^p = V-u innx(0, T), (1.8) 

with the appropriate boundary and initial conditions. 

A priori error estimates for solving (ll.7|) - p.8|) were obtained in [510 El [TH] . In m, Geveci derived 
L°°-in-time, L^-in-space error bounds for the continuous-in-time mixed finite element approximations 
of velocity and stress. In [610, a priori error estimates were obtained for the mixed finite element 
approximation of displacement which requires less regularity than was needed in [9]. Stability for a 
family of discrete-in-time schemes was also demontratred. In m. an alternative mixed finite element 
displacement formulation was proposed reducing requirement on the regularity on the displacement 
variable. For the explicit discrete-in-time problem, stability results were established and error esti¬ 
mates were obtained. The effectiveness of the method analyzed in m was demonstrated in [12] by 
providing simulations using both lowest-order and next-to-lowest-order Raviart-Thomas elements on 
rectangles [25] . 

The purpose of this paper is to analyze an implicit time-stepping method combined with the mixed 
finite element discretization proposed in [13]. We prove the stability of the proposed method by using 
energy estimation, and show in particular that it conserves certain energy. We also invertigate the 
convergence of the method and prove optimal a priori L°^{Lp) error estimates for both displacement 
and pressure. The rest of the paper is organized as follows. In sections 2 and 3, we introduce notations 
and describe the weak formulation of the problem. The fully discrete mixed finite element method is 
presented in section 4. Stability results are established in section 4 and optimal a priori error estimates 
are obtained in section 5. Conclusions are given in the last section. 

2 Notation 

We shall use the following inner products and norms in this paper. The L^-inner product over 14 is 
defined by 

(u, v) = / uv dfl, 

Jn 

inducing the L^-norm over 14, = (u,u)^/^. The inner product over the boundary 914 is 

denoted by 

{u, v)= uv dI4 

Jan 
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for u, V € ff2+^(r2) with e > 0. We introduce the time-space norm: 

lbllL2(0,T;L2(n)) = ll^'llL2(i2) = 11^1 . 

The time-space norm || • \ \ioo(^i 2 'f is similarly defined. In addition to the spaces, we use the standard 
Sobolev space for mixed methods: 

H{n, div) = {v:ve V ve 

with associated norm 

ll^ll//(o,div) = ll^llL2(n) + l|V •'w||L 2 (f^), 

where 

( m 

Y.HWhin) 

i=i 

For the time discretization, we adopt the following notation. Let be a positive integer. At = 
T/N, and = nAt. For any function v of time, let u” denote u(t”). We shall use this notation for 
functions defined for all times as well as those defined only at discrete times. Set 

^ ’ 

^n-,e ^ Q^n+1 ^ ^ _ 20)yn ^ Q^n-1 ^ 

where 0 < 0 < 1. We also define the discrete /°°-norm for time-discrete functions by 



l'*^IU?>,(0,T;L2(n)) — ||u||;oo(i2) 


max m 
0<n<N 


‘llL2(n). 


3 Weak Formulation 


The finite element approximation of the wave problem is based on its weak formulation which is 
derived in the usual manner. Integrating by parts and using the data on the boundary of fl, we obtain 
the weak formulation m- For any t > 0, find {u{t),p{t)) G V x IT such that 


{u{0),v) 

= (u^,v) VveV, 

(3.1) 

(«t(o),-y) 

= (v^,v) VvgV, 

(3.2) 

(A“^p(0),u;) 

= (V-u°,w) VweW, 

(3.3) 

{puu{t),v) + {p{t),V -v) 

= {f{t),v) Vw G T, t > 0, 

(3.4) 

{X~^p{t),w) - {\/ ■u{t),w) 

= 0 Vw G IT, t > 0, 

(3.5) 


where V and IT are given by 


V = {v € H{Q, div) : • i/|rjv = 0}, 
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W = H 2'''^(r2) for any e > 0. 

The present formulation requires less regularity on the displacement than standard approaches. For 
instance in mm it is necessary that Vp G fr(0,div) so that V •« G Here, it is only required 

that V ■ u & H 2 ^ and it can be verified that the solution u of problem (ll.lD - (ll.5p with p = AV • u is a 
solution to dll-dSS]) , see m- 

Differentiate ()3.5p with respect to time to obtain 

{X~^Pt,w) — {V ■ Ut,w) = 0 Vw^W. (3.6) 

We next assume / = 0 and choose v = Ut and w = p in ()3.4p and ()3.6p . respectively, so that 

{puu,ut) + {p,\/-Ut) = 0, 

{X~^pt,p) - {V ■ut,p) = 0. 

By adding the two equations, we hnd that 

{pUtt,Ut) + {X~^Pt,p) = 0, (3.9) 


(3.7) 

(3.8) 


or 


2 dt 


p2Ut 


2 1 d 

L2(r2) 2 dt 


X 2p 


L2{Q) 


= 0 . 


Thus, in the absence of forcing, the (continuous) energy 


1 

2 


1 


p2Ut 


2 1 

+ TT 

L2(o) 2 


_ 1 

A 2p 


2 

L2(o) 


(3.10) 


is conserved for all time. It will be shown that a similar form of energy is conserved by the numerical 
solution of the wave problem. 


4 Finite Element Approximation 

For the finite element approximation, we let {Sh}h>o be a quasi-uniform family of finite element 
partitions of D, where h is the maximum element diameter. Let 17/i x Wh be any of the usual mixed 
finite element approximating subspaces of V x W, that is, the Raviart-Thomas-Nedelec spaces [T91I23] . 
Brezzi-Douglas-Marini spaces [1] , or Brezzi-Douglas-Fortin-Marini spaces [3] . For each of these mixed 
spaces there is a projection H/i : H{Xl,div) -^Vh such that for any 2 : G H{Q,div) 

{X/■ UhZ, w) = {V ■ z, w) yw£Wh- (4.1) 

We have the property that, if z G div) n/f^(D), then 

WUhZ - z\\o <Chmz\\j, l<j<k, (4.2) 

where k is associated with the degree of polynomial and || • ||s is the standard Sobolev norm on 
(iL*(D))™. Here and in what follows, C is a generic positive constant which is independent of h and 
At. 

For (j) G W, we denote by Vhff) the L^-projection of (/> onto Wh dehned by requiring that 


{Vh(l), w) = {(j), w) Vrc G W/i. 


( 4 . 3 ) 
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If (/> G VF n then we also have 

\\Vh(l>-<P\\s<Ch^-^\<P\\j, 0<s<k, 0<j<k. (4.4) 

The semidiscrete mixed finite element approximation to {u{t),p{t)) is to seek {U(t), P{t)) G Vh x Wh 
satisfying 


{U{0),v) 

= (n/in°,r;) 

Vw G Vh, 


(4.5) 

ium,v) 

= {IlhV°,v) 

Vv G Vh, 


(4.6) 

{P{0),w) 

II 

(o 

yw G Wh, 


(4.7) 

ipUtt{t),v) + {P{t),V-v) 

= ifit),v) 

yv G Vh, 

t > 0, 

(4.8) 

{X~^P{t),w) - {\/ ■U{t),w) 

= 0 

yw G Wh, 

t > 0. 

(4.9) 


Existence and uniqueness of a solution {U{t), P{t)) to the variational problem (|4.5I) - (|4.9D is shown in 

[13]. 

The fully discrete mixed finite element 0-scheme is then defined by finding a sequence of pairs 
{U^,P^) eVhxWh,0<n<N, such that 


{U^,v) 

(P^w) 


pdtU 2,w) + (dtP^,V v 


+ ^(P»,V.v) 


{pduU^,v) + iP^-^^,V -v) 
^^-lpn+1/2^ U;) - (V • 

Equation (I4.12p is derived from the following expansion: 


\/veVh, 

{p°,w) Vw€Wh, 

(^^f + eAt^BjKv^ 

\/v€Vh, 

(r’^u) yvGVh, 

0 \/w G Wh- 


=u^ P Atv^ + Ar 


+ 


-0U' 


+ 0(At^ 


(4.10) 

(4.11) 


(4.12) 

(4.13) 

(4.14) 


The present 0-scheme is explicit in time if 0 = 0 and implicit otherwise. The existence and uniqueness 
of a solution to the resulting linear system for a nonzero value of 0 follows from the unisolvancy of the 
mixed formulation of the following elliptic problem: 


V • (AV0) + 


0Ab 


p(/> 


(p 


0 in n, 

0 on 


The explicit case has been considered in [T3|. As expected from an explicit scheme, the method is 
conditionally stable. As a stability constraint, it requires to choose 


At = 0{h). 


In the next sections, stability and convergence properties of the proposed 0-scheme are analyzed. 
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5 Stability Analysis 

We derive sharp stability bounds based on the energy technique and show that the proposed scheme 
conserves certain energy. We consider (j4.13p and (I4.14p for the homogeneous case 


+ = 0 yveVh, 

^^.ipn+i/2^w) - (V • w) = 0 Vw€ Wh- 


(5.1) 

(5.2) 


We will make use of the inverse assumption, which states that there exists a constant Cq independent 
of h, such that 

l|V • <)>||l2(o) < Coh~^\\4>\\L'2(n) 
for all (j) € Wh- The following stability result holds. 


(5.3) 


Theorem 1 . The fully discrete scheme (I4.inh - (l4.14p is stable if 




and conserves the discrete energy 


-2 




-^ ) ||A-t5tP"+2||2 + 


- 2 pn-\- 2 I 


(5.4) 


(5.5) 


pn\9 

'2At 


The scheme is unconditionally stable if 0 > 1/4. 

Proof. If we subtract (j5.2p from itself, with n + 1/2 replaced by n — 1/2, we find that 

(y-l(pn+l _ _ (y • {U^+^ - U^-^),w) = 0. 

As (15.ip holds for all v € V/i and (15.6p holds for all w € Wh, we choose v = dtU"^ and w = 

(pa^^^7^9^c/") + (p”-^v•a^c/") = 0, 

■dtU^,P^'^) = 0. 

By adding (15.7p and (15.8D we obtain 

{pdttU\dtU^) + {x-^dtP^,p^'^) = 0. 

Note that 

pnfi ^ /Xt^eBuP"^ + P'^ 

Hence, (15.9p can be rewritten as 

{pdttU^,dtU^) + M^ fe- {X-^duP'^,dtP^) + \ {x-\P^+"2+p--i)^BtP^] = 0. (5.11) 


(5.6) 
so that 

(5.7) 

(5.8) 


(5.9) 


= /xe(e--\ BttP'^ + - (+ P"“^ 


(5.10) 


Using that 


B un ^ ^ - Q ^n ^ ^ 


At 
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we deduce that 


{pBttU^,BtU^) = ^ {pBtU^+-2 -pBtU^-KBtU^+'^ + BtU^--^) 


1 


2At 


{pBtU'-^+KBtU^^-^) - {pBtU^-^BtU^-'^) 


and similarly 


i-ia r>n a D^^^ _ ^ T/x-iu. 0^1+^ 


{X-^BuP^,BtPn = ^ [{\-^BtP^+i,BtP^n) - {\-^BtP^-"2,BtP^-'2) 


We also have 


a-1(p«+ 2+p^-t),atP") = ^(A-^p"+t+ A-^p’"-t,p^+t 


^ lon+i I \—IdJ^—^ TDn+h 


At 
1 

At L 


(A-lpn+i^pn+i) _ 


Hence, ()5.1ip is equivalent to 


^ (<*’ - <■’ I =». 


where Ej^ ^ is the quantity defined by (15.5p . This relation indicates that Ef^ ^ is conserved for all 

n+t 

time, which guarantees the stability of the scheme if and only if E^^ ^ defines a positive energy. A 
sufficient condition is that 




+ AP e- 


1 


A-5 9iP"+i 


> 0 


for all n > 0. Clearly, the scheme is unconditionally stable when 9 > 1 /4. Now, using Cauchy-Schwarz 
inequality and the inverse assumption (|5.3I) . we obtain 


X-^dtP^+^,w] = {V-dtU'^+^,w 

1 


< 


< 


V • ac/”+t 


Co 


BfU^+2 


L2(n) 


L2(n) 


\w 


llL2(n) 

llL2(n)- 


By setting w = dtP'^^^, we see that 


x-^dtP'^^^ 

2 ^ ^ 

BtU'^+^ 


dtP^+h 


L2(f2) “ h 


L2(n) 



or 


< 


X-^tP^+l 


CoXl 


hpt 




L2(o) 


< 


hpl 


t r\ yTTI-\- — 

p^dtU ^ 


L^{n) 

x~^BtP^~^^ 

L2(r2) ■ 


L2(0) 


Hence, a sufficient condition for stability is given by 


\\pWtU'^+^\f + At^ [9- 


Oln 


\\pWtU^+2\\^ > 0 


4y h?Po 


which completes the proof. 


□ 


The case with 0 = 1/4 is interesting because the form of the discrete energy in this case is similar to 
that of the continuous problem. In addition, one can verify that the time truncation error is minimized 
over the set of all 0 > 1/4 when 9 = 1/4. 
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6 Convergence Analysis 


In this section, we prove optimal convergence of the fully discrete finite element solution in the L°°(L^) 
norm. Some of the techniques used in the proofs can be found in previous works [lade]. In order to 
estimate the errors in the finite element approximation, we define the auxiliary functions 


= U^- UhU^, ri^ =u^- UhU^- 


^pn_ 


C = p"- w, 


where and Vh are defined in Section 4. From (|3.4p - (|3.5|) and (|4.13I) - (I4.14|) . and the properties of 
the projections II/i and Vh, we arrive at 

{pdttX'^,v) + {C"'^,V-v) = {pdu‘n'^,v) + ,v) 'iv h, n > 1, (6.1) 

(^_l^„+l/ 2 ^^) _ (V.^n+i/ 2 ^^) ^ ^weWh, n>0, (6.2) 

where r"' = p{u'^f — BuU^)- Another equation has to be derived for the initial errors and 
Consider (13.4p at n = 0 and n = 1, respectivey, and subtract the resulting equations so that 


A use of Taylor’s formula with integral remainder yields 


pdtuf^,v + {dtp 2 ,V-v) = {dtf^,v 


,d^u 


o — n 0 

dtU2 =v^ + —< + 


1 


r-Ai 


2AtJ, 


(6.3) 


(6.4) 


Using (16.3p and p6.4K . we readily obtain 


pdtU 2 ,v) + 9At‘^ (dtp^ ,V ■ v) = 9At^ (dtf^ ,v) — OAt'^ [ pdtu^^,v 


+ {pv°,v) + ^{pu^t,v) 


2At Jo 

from (|4.12l) and taking into account (|4.1ip and 
At 


Ai / 

+ TT:rT / (At - tf ( Pp-^,v ] dt. 


1 


(6.5) 


At / q3,, 

PlVT,v]dt. 


( 6 . 6 ) 


dt^ 

Subtracting (16.51) from (14.121) and taking into account (14.111) and (13.41) to arrive that 

^^pdtxKv^ +9At^ -v^ = (^pdt‘n^,v^ + {p{UhV^ - v°),v) 

" ‘ " " ^— / l /\i — i\^ \ n 

dt^ 

Note that = 0 and = 0. We now state and prove our convergence result. 

Theorem 2. If u G L°°{H{n-, div)), ^ € L^{L^{n)), ^ € L°^{L^{Vt)), and p € then 

for {!/", P"} defined by (l4.1UI) - p4.14p there exists a constant C independent of h and At such that if 


A,2,1 


(6.7) 


then the following a priori error estimate holds: 


p2(u — U) 




+ 


A ^{p — P) 

< C{h^ + At2) 




d^u 




L°°(L2) 


+ \ \p\\l^{l^) 2 


( 6 . 8 ) 


where r is associated with the degree of the finite element polynomial. 
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Proof. We first rearrange (|6.1D in the form 

{pduX^,v) + Af (e - {duC, V • t;) + i V • v) = ipBuri^, v) + {r^,v). (6.9) 

Summing over time levels and multiplying through by At yields 


ipdtX"'^^^ - pdtX^v) + At'^ f f9fC+2 -v 


+ 


2,V -u) = (^pdtT]'^+2 -pdtvKv) + ( At^r*,' 


2=1 


2 = 1 


Upon defining 


we verify that 


n—1 


= 0, 0- = At^r+i 


i=0 


At 1 At 




2 = 1 


Taking into account (16.6p and that dtf'^ = , (|6.1Up becomes 


where 


.R” = At ^ r* + p{UhV^ - v^) + OAfpdtU 


f-At 


2=1 


2At 


(Ai - 0 


( 6 . 10 ) 


(p9iX"+^,t») + At2 (^0 - 1) V • v) + ■ v) = (pdtV^^Kv) + iR^,v), (6.11) 


Since = ^"'’^ 2 , (16.2p reads 


x-^dt(tP+^^,w] - h-x'^^Kw) = (x-^c^^,w 


( 6 . 12 ) 

Choosing v = and w = in po.lll) and (I6.12p . respectively, and adding the resulting 

equations, we arrive at 


{pdtx"+Kx"^-^)+At^ y V •x"+^) + [x-^dtr^Kr ^'2 

—l^n+i ^n+i 




Again, choose w = in p6.2p so that 




Substitution into (|6.13p yields 


(/o9tx”+ix”+^)+At" ( ^ - 7 ) (A-i5ir+t,r+^) + {x-^BtB"^i 






=AP 


(A"'r+i + (p5ir7"+Ux"+5) 




(6.13) 


(6.14) 
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The terms on the right-hand side of (j6.14li are bounded using Cauchy-Schwarz inequality as 

1 


A-l^n+i^^^^n+f 

< 

A-ic^+i 

L2(n) 


< 

pdtrf^^ 

L2(o) 


< 1 

lL2(n) 

^n+i 


< 


L2(f2) 


\-2dte^2 
L 2 (o) 


L 2 (r 2 ) 




L 2 (n) 

A-4"+^ 


L 2 (n) 


We now distinguish the cases where 0 > 5 and 9 < ^. In the first case, we sum on (|6.14|1 over time 
levels and multiply through by 2 Af. This results in 

2 

L 2 (n) 



2 

- 0 

2 

-|- 

A-i</>"+i 

2 

A 2())0 


L 2 (n) 


L 2 (Q) 


L 2 (o) 



+[ 9 - 


1 




< 2 Af^ 9 - 


I El|A-h- 


i=0 


+ 2Af^||x'+5 


i=0 

n 


L 2 ( 0 ) 


L 2 (n) 


L 2 (n) 


pdtrf^^^ 


A-5^° 


L 2 ( 0 ) 


— T ti-l-1 


A 2 ^ 


-h 2At^ A"5(^*+5 


i=0 



A-if+i 

L^(n) 



L 2 (o) 

L2(o) 


L2(o) 


+ 


A-if 


L 2 (o) 


(6.15) 


L2(n) 


Since 

A-if 

< 

A 2 ^ 

and 

/ 02 X 2 

< 




L^{n) 




L 2 (n) 



Z“(L 2 ) 


, then 



2 

-l- 

A-i())'^+i 

TAf/'e-i') 

A-ir+^ 


L 2 (n) 


L 2 (o) 4y 



2 

L 2 (o) 


< 4Af 9 - 


+ 


2At 


/’o 


+ 2M 


p'^X 


A 2 


A 2 ^ 




i=0 


/ n 


li 




vi=0 
n 


L 2 (r 2 ) 


L 2 (n) / 

n 

+ EII^‘ 


(6.16) 


j=0 


Il 2 (o) 




^ A 2 f +2 


^*=0 


L 2 (f 2 ) 


Applying the algebraic inequality: ab < |a^ + to the right-hand side of (|6.16p shows that 



2 

-|- 

A-i())^+^ 


A-i^”+i 


L2(n) 


L2(C!) V 47 



2 

L 2 (o) 


1 


<-Af 0 - - 


1 


A 2 ^ 


Z“(L 2 ) 


Af-1 


+ 8(0--) At ^ A-if+^ 


i=0 


L 2 (n) 


1 

+ 2 


1 

+ 2 


/N-l 


N-1 


P^X 


A 21 ^ 


Z°°(L 2 ) 

2 


Z°°(L 2 ) 


/ AT-l 

-h4 ( At ^ llA^if'^i 


(6.17) 


i=0 

2 


lL 2 (!rZ) 


i=0 


L 2 (o) 
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If we take the supremum over n on the left-hand side and use the fact that NlS.t = T, we conclude 
that 

2 


1 

2 

1 

2 

1 


+ 

A 2 

< (7 

A- 2 ^ 






/N-l 

+ C/S.t^ I ^ I 
. 2=0 


L 2 (n) 


/N-l 

+ CAt^ X] 11 ^* 


(6.18) 


. 2 = 0 


iL 2 (Q) 


For the case 0 < we can follow the analysis presented in [laile] to derive error estimates similar 
to ()6.18l] under condition (j6.7p . 

To complete the proof, we need to bound each term on the right-hand side of (I6.18p . The first 
term can be bounded using the approximation properties. Similarly, we have 

Af-l 




2 = 0 


d^u 


dt^ 


Ll(0,T;L2(n))^ 


For the last term on the right-hand side of ()6.18p . we have 
N-l 

AtY^m^^n) < C\\R\\i^^l^) 

N-l 

CAtYn^m + cMu^v^ 

i=l 

+ceAt‘^ 


2=0 


< 


1 


1 n9 <9^1* 

P5t«"t 

+ c 

L2(r2) 

2AfX ') atW-" 


L 2 (o) 


To estimate we make use of the identity 

1 


= < + 


6Af^ 7_At 


■At q4,. 


(6.19) 


From the Taylor’s expansions of and ^ about u'l’t', 


/■At o4, 

<+' = u’i + AtuS, + / (At - |4)7=-(t” + ») ds, 


Jo 


and 


we obtain 


u 


dt^ 

d^u 


/■O //4» 

= < - Atul, + j JAt - \s\)-^{r + s) ds, 


/ /\ -f- 

JAt-\a\)-^{e + s)ds. 


^tt ~ '“tt 


( 6 . 20 ) 


Subtracting (I6.19P from (16.201) yields 

. 1 
u'tt 


- Buu^ = 


GAr /_ 


/■At o4„ /-At q4, 

LJ‘' - V(‘" 


ds. 


Hence, 


|2 

lL 2 (n) 


= I|P 


n:9 


U, 


- dttvJ 


|L 2 ^ 


< CAr 


f At 


'-At 




+ s) 


ds < CAt'^ 


L2{n) 


1 cBu 


L°°(L 2 ) 
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and therefore 


AtY,nHn)<CAt^ 


2=1 

Similarly, we have 


1 d‘^u 


V At < CTAt^ 

L°°(L 2 )^ 


d^u 




pdtult 


and 


2 

L 2 (o) 

r-At 


A(/ 


< CAt^ 


L 2 (o) 


/•At 

i d^u 

L 



(it < CAt^ 


L 2 (r 2 ) 


9t3 


L°°(L 2 ) 


1 \ 


dt 


< CAr 


i-At 

1 d^u 

L 



dt < CAt^ 


L 2 (n) 


1 ci^u 


L 2 (n) 

Finally, using the approximation property (|4.2p and combining all the bounds, we arrive at 

TV-l 


L°°(L 2 ) 


^ ll■^*llL2(o) < Cp/i^ + At^), 


i=0 


which completes the proof of the desired estimate. 


□ 


Remarks. It is worthwhile to mention that the time discretization method is fourth-order accurate 
when 0 = 1/12. To preserve the temporal accuracy of the hnite element scheme one has to modify 
(14.121) carefully to obtain an appropriate initial value U^. The analysis presented in [15] can be used 
to derive optimal a priori error estimates in this case. 


7 Conclusions 

We proposed and analyzed a family of fully discrete mixed finite element schemes for solving the 
acoustic wave equation. We derived stability conditions for conditionally implicit stable schemes 
covering the explicit case treated by Jenkins, Riviere and Wheeler m- The error estimates established 
provided optimal convergence rates for the use of mixed finite elements methods in solving the acoustic 
wave equation. 
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